Streaks to Rings to Vortex Grids: Generic Patterns in Transient Convective Spin-Up 



J.-Q. Zhong, 1 M. D. Patterson 1 ' 4 J. S. Wettlaufer 1 ' 2 ' 3 
1 Department of Geology & Geophysics, 2 Department of Physics and 
3 Program in Applied Mathematics, Yale University, New Haven, CT 06520-8109, USA and 
4 Faculty of Engineering & Design, University of Bath, Bath, BA2 7 AY, UK 

(Dated: June 29, 2010) 

We observe the transient formation of a ringed pattern state during spin-up of an evaporating 
fluid on a time scale of order a few Ekman spin-up times. The ringed state is probed using infrared 
thermometry and particle image velocimetry and it is demonstrated to be a consequence of the 
transient balance between Coriolis and viscous forces which dominate inertia, each of which are 
extracted from the measured velocity field. The breakdown of the ringed state is quantified in terms 
of the antiphasing of these force components which drives a Kelvin-Helmholtz instability and we 
show that the resulting vortex grid spacing scales with the ring wavelength. This is the fundamental 
route to quasi-two dimensional turbulent vortex flow and thus may have implications in astrophysics 
and geophysics wherein rotating convection is ubiquitous. 
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The impulsive rotation, or change in angular velocity 
of bounded fluid described by the Navier- Stokes equa- 
tions results in the propagation of stresses into the inte- 
rior. Despite the ubiquity of the process, it remains suf- 
ficiently complex to insure its centrality as a long stand- 
ing and fundamental problem in fluid mechanics [TJ |2]. 
Fluids are often heated or cooled at their bounding sur- 
faces and whether they are homogeneous or stratified, 
there are common conundrums that influence a range 
of disciplines from astrogeophysical flows to engineer- 
ing phenomena [e.g., 3-8 . Here, we examine pattern 
formation during transient spin-up experiments where 
buoyancy forcing is supplied by evaporative cooling from 
the top surface of a layer of water of depth H that is 
probed using infrared thermometry and particle imaging 
velocimetry (PIV). It is found that this localized sur- 
face cooling drives convection that evolves into a radi- 
ally symmetric ringed, or "bulls-eye", pattern near the 
free surface (Fig. [TJ}) itself decaying vertically creating 
a "bowl" shaped region (Fig. [5J. We extract velocity 
fields to uncover the origin of the transient, symmetric 
and depth-dependent ring pattern. It is demonstrated 
that this results from a temporary balance between the 
Coriolis and viscous forces which dominate inertia during 
an Ekman spin up time te = ^H 2 /vQ of order unity, 
where v is the kinematic viscosity. Finally, we find that 
beyond this transient state, the ring pattern breaks down 
to the canonical quasi-two dimensional vortex grid state 
[9l [10] in a manner reminiscent of a Kelvin-Helmholtz in- 
stability (Fig. [lj). It is found that the number density 
of vortex cores TV is inversely proportional to the square 
of the ring wavelength A, itself decreasing with Q. 

The patterns that emerge in rotating convection dis- 
play the competition between Coriolis and viscous forces 
across a range of space and time scales [e.g., [5j [8]. A sig- 
nature process in rotating convection contrasts the im- 
plications of Taylor-Proudman theorem, which forbids 
vertical motion for an inviscid constant density rotating 
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FIG. 1: A sequence of particle-number-density images show- 
ing the flow field at a horizontal level 0.6 cm below the free 
surface of a volume of water H — 11 A cm deep with a 22.9 
x 22.9 cm cross section. Unless otherwise specified dimen- 
sionless time is t = t' /te, with time — t! . (a) t = just 
before the initiation of rotation, (b) t — 2.6, the flow field 
was axisymmetric consisting of regular concentric rings, (c) 
t — 3.9, wavy billows appeared upon the breakdown of the 
rings, (d) t = 5.0 the vortex grid state. In this experiment 
Ro = 1.4xl0 -3 and r E = 167 s. The movie is online 20 . 



fluid, with the buoyancy forcing, which can drive narrow 
regions of weak convection in fluids with finite viscosity. 
Thus, the formation of vortices or sheets in the interior 
of a rotating fluid subject to buoyancy forcing can be 
characterized by the flux Rayleigh number Ra describing 
the strength of buoyancy transfer in the fluid, the Rossby 
number Ro capturing the relative influence of buoyancy 
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FIG. 2: Thermal images of the surface taken at four stages 
during the spin- up process. In the experiment Ro = 5xl0 -4 , 
the rotation period is 5.5 s and hence te = 119 s. The images 
were taken at (a) t = 0, (b) t = 2.9, (c) t = 3.7, (d) t = 8.7. 
The color bar is the temperature scale in °C. 

forcing to rotation and the Taylor number Ta character- 
izing the relative influence of rotation to viscous forces: 

BH 4 n r~B~ m /2QH 2 \ 2 
Ra =^7> R °=yWlP' Ta= {—) ' (1) 

where B is the buoyancy flux and k is the thermal dif- 
fusivity of the fluid (in our case water). It is well known 
that in rotating and non-rotating convection [e.g., [8] the 
critical Rayleigh number (Ra c ) for the onset of convec- 
tion depends on the nature of the thermal boundary con- 
ditions (fixed temperature, fixed flux, or mixed) and that 
in the former case whether both bounding surfaces are 
free, rigid or one is free and one is rigid; Ra c converges to 
a large Taylor number scaling of Ra c — » constant Ta 2 / 3 
[8]. Far beyond critical, a particularly compelling ex- 
ample of convective pattern formation is the transient 
ring state observed during spin-up. The present study, 
in combination with previous work [9j [10] , confirms that 
this transient state, and its breakdown into a vortex grid 
is generic in that it occurs under all thermal and me- 
chanical boundary conditions described above, as well as 
in both cylindrical and square cross section cells. Our 
method quantifies the underlying causes for the evolu- 
tion of the entire process. 

The square cross-section cell containing deionized wa- 
ter initially at rest is set into rotation about its ver- 
tical central axis to values of from 0.1 to 1.0 rad 
s _1 . The bottom and side walls are insulated and main- 
tained at the same temperature T w = 30 ± 0.1° C and 
the relative humidity above the fluid is held constant at 
35% ± 1% during the course of an experiment. Thus, 
evaporation induced a constant heat flux Q = 56.0 W 



m 2 driving a buoyancy flux B — 3.7 x 10 8 m 2 s 3 , so 
that Ra = 3.3 x 10 8 > Ra c , 5 x 10~ 4 < Ro < 10~ 2 and 
2 x 10 7 < Ta < 1.4 x 10 9 , thereby insuring a rotationally 
dominated regime of flow. 

The number density of seed particles in the flow is cap- 
tured using a grid-based PIV system [11 to image the 
flow structure as a function of time. The velocity data 
were calculated in a 1000 xlOOO grid with a spatial res- 
olution of 0.2 mm/pixel. For Q = 0, Fig. [TJa) shows 
a random network of thin sheet-like downwelling cold 
streamers typical of evaporatively-induced turbulent con- 
vection [12]. Shortly after rotation begins, these down- 
welling plume structures are sharply modified and orga- 
nized into concentric rings. The rings form first at the 
center, then sequentially grow to cover most of the hori- 
zontal plane by t = 2.6 rotating counterclockwise with 
uniform azimuthal velocity in a quasi-stationary state 
(FigjlJ)). Kelvin- Helmholtz billows appear on the rings 
at t = 3.9, grow and roll up into counterclockwise vortices 
(FigjlJ). Finally, at t = 5.0, the azimuthal symmetry of 
these vortices is lost as they assemble into the regular 
vortex grid seen in Fig[TJd). At this stage we imaged the 
vertical fields to find classical rotating convective colum- 
nar structure plfTSUTTj. 

Infrared thermometry confirms the hypothesis that the 
basic dynamics of ring formation and breakdown is the 
same with free surface evaporative cooling and bottom 
heating [9], Rayleigh-Benard [10 boundary conditions, 
and free surface evaporative cooling with insulated bot- 
tom boundary conditions as studied here. The surface 
temperature maps (Fig. [2| reveal the full evolution of 
the thermal forcing and associated pattern formation, 
and a positive correlation between thermal fluctuations 
and vertical velocity p~8j [19] . The evolution of the pat- 
terns seen in Figj2|a-d) are to be compared with those 
in Figjlja-d). Although these are two separate runs, the 
evolution is essentially the same; the overall temperature 
drop from the wall to the center is AT w 2.0°C and is 
the same at different stages of spin-up. The width of the 
vertical thermal boundary layer at the wall is St ~ 10 
mm, evolving slowly in time. 

Fig.Qb is the map for the maximal ring structure hav- 
ing a radial temperature amplitude of oscillation ST = 
0.2°C = 0.1 AT and a minimal St- The radial oscillation 
between descending (dark) and ascending (light) fluid is 
clearly seen, finally merging with the thermal boundary 
layers at the walls. When the Kelvin-Helmholtz insta- 
bility is initiated in Fig.Qc the ringed-shape thermal 
plumes break down into vortices, the centers of which 
contain sinking fluid. 

Because the flow is rotationally dominated we complete 
our analysis by considering the azimuthal component of 
the momentum equation in a co-rotating reference frame 
in cylindrical-polar coordinates, 

^r + [(«-v)tf| fl = [vv 2 u] 9 -2[nxu\ 9 , ( 2 ) 
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FIG. 3: Fluid accelerations measured in the experiment with 
Ro — 1.4 xlO -3 . (a) Averaged accelerations as function of 
time. Black curve: blue curve: (|/ c |); and red curve: 

The vertical dashed lines delineate when the ringed 
structure is maximal, t — t\ — 2.6, and when the rings started 
to break down t — £2 = 3.1. The inset focuses on £1 < £ < £2- 
(b) Fluid accelerations as functions of the radial distance at 
three stages. Blue curve: f c ; red curve: f u . (c) Correlation 
function C(t) between f c and f u as a function of time. Vertical 
bars denote fluctuation amplitudes. 



where u = (u r , uq) describes the radial and azimuthal ve- 



locity components and 



A 



denotes the azimuthal com- 



ponent of A. We nondimensionalize using the Ekman 
length Le = yjvj^l, and te to extract the Coriolis (/ c ), 
viscous (f u ) and inert ial (fi) forces per unit mass from 
the PIV velocity fields by integration of each term of Eq. 
Q over = (0, 27r). The areal averaged forces per unit 
mass for a component fj are defined as 



<i/,i> 



\fj\rdr, 



(3) 



from the ring center (r = 0) to the maximum radius 
(r = r m ). 

The time evolution of (|/ c |>, (\f„\) and <|^|> (Fig.(§a) 
shows that upon spin- up the fluid inertia (\fi\) decreases 
most rapidly and (|/ c |) converges with (\f u \) as the ring 
state begins to appear at the center. Most striking is the 
abrupt appearance of the ring pattern which fills the cell 
surface at t = t\ and then breaks down at t = £2 (vertical 
dashed lines). Hence, the dissection of the velocity field 
shows that the origin of the transient ring state is that 
during h < t < t 2 , (\f c \) = (\fu\) - 10(|/i|>; the Coriolis 
force balances the viscous force both of which dominate 
inertia. 




FIG. 4: Flow structures revealed by the thermal images 
for different Ro. From (a) to (f) Ro=1.2x 10" 2 , 4.0 xlO -3 , 
2.2xl0 -3 , 1.4xl0" 3 , 7.6xl0 -4 and 5.0xl0 -4 respectively. 
With the exception of (a) the image was taken when the 
ringed structure was maximal. Image (b) is taken at the 
threshold value Ro c for ring formation. 



Figs.([3|b,c display the radial dependence of f c and f v 
at three times and their correlation function 



C(t) 



f-J m {fc-{fc)){h-(h))rdr 



yJUZ r o m (fc ~ {fcWrdr] [fZ m Uv ~ </„» 2 rdr] 

(4) 

as a function of time, where the (fj) are defined by Eq. 
([3| but without taking the absolute value. Thus, the 
maximal ring pattern of Fig.Qb appears as a maximum 
in C(t) at t = ti when f c = f v and they are in phase as 
seen in the middle panel of Fig.([3|b. That the ring pat- 
tern breaks down at t = £2 through a Kelvin-Helmholtz 
instability is demonstrated by the abrupt drop in C(t) 
signaling the loss of phase between f c and f v at t = 3.1 
and the beginning of rotational dominance in the hori- 
zontal force balance driving shear and vortex formation. 

The temperature maps in Fig. Q show the flow pat- 
tern at different Rossby numbers when the ringed struc- 
ture is maximal. The ring wavelength A decreases with 
increasing Q and there is a clear threshold Rossby num- 
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FIG. 5: Particle- number-density images versus fluid depth z 
in cm. (a) z — 0.6, (b) z — 3.1, (c) z — 5.7 and (d) z — 8.3. 
Here Ro = 1.4xl0 -3 the images are taken at t = t2 when 
ring breakdown begins. 

ber, here Ro c = 4.0xl0 -3 , above which the ring patterns 
are absent (Fig. |4|i). Indeed, for small Q horizontal mo- 
tion is too weak to suppress vertical thermal convection; 
the flow resembles the ft = case of Fig. ((2j)a- When 
t^>t2 the ringed pattern evolves into a regular vortex grid. 
After spin up we systematically varied Vt and measured 
the number of vortices per unit area A/", and the mean A 
at t = t\ . The average spacing between the rings (FigfTj}) 
is related to the mean distance between adjacent vortex 
cores (FigJTJi) according to A = C/>/~N, with a best fit of 
C = 0.605 ±0.015. 

Finally, to demonstrate the vertical extent of the 
ringed state we present the flow pattern simultaneously 
at different horizontal levels at the t = £2 just when the 
surface pattern begins to break down (Fig. |5|. Here we 
see the "bowl" structure of the ringed pattern with the 
central ring disappearing beyond a depth of z = \H be- 
low the surface (Fig. [5]i). Because of the finite lifetime of 
the rings, forming at the center and moving radially out- 
ward, the depth of fluid penetration necessarily decreases 
with radius. 

In summary we have used a combination of PIV and 
infrared thermometry to quantify the evolution of an 
evaporatively driven convecting spin-up system. We find 
that, beginning at 2.6 t#, a balance of Coriolis and vis- 
cous forces dominating over inertia leads to a transient 
concentric ringed convection pattern with radial sheets 
of downwelling and upwelling fluid. When these domi- 
nant forces lose phase coherence, the sheets breakdown 
into the vortex grid typical of quasi-geostrophic turbu- 



lence. Moreover, the ring spacing and the vortex core 
spacing after ring breakdown obey a simple scaling re- 
lation. Our methodology demonstrates the generality of 
the ringed pattern state. As a practical matter, the great 
sensitivity found here of convection driven by evapora- 
tion under room temperature conditions suggests that 
many experiments with an exposed upper surface will 
display unintended dynamical effects. Finally, we deter- 
mined the critical Rossby number for the existence of a 
ringed state and this suggests possible implications for 
both astrophysical and geophysical convection wherein 
rotating convection is ubiquitous. 
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